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Abstract 



We consider a mass-conserving bistable equation with a saturating flux on an interval. This 
is the quasilinear analogue of the Rubinstein-Steinberg equation, suitable for description of 
order parameter conserving solid-solid phase transitions in the case of large spatial gradients 
in the order parameter. We discuss stationary solutions and investigate the change in 
bifurcation diagrams as the mass constraint and the length of the interval are varied. 
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1 Introduction 

Consider a bounded interval, f2 = (0, L) C M, L > 0. To extend the Ginzburg-Landau solid-solid 
pliase transition theory to high gradients in the order parameter u, Rosenau [7] proposed the 
following free energy functional: 



where the diffusion coefficient e G (0,oo), the interface energy "^{s) is a convex function of its 
variable that grows linearly in s so that, for example. 
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and W{u) is the bulk energy which we take to be the double-well function, 



Win) = ^-^. 




(1.2) 



is the well-known mean curvature operator and can be viewed as the flux which satisfies 



The bistable Rosenau equation (11. 2p with the physically relevant Neumann boundary conditions 
is studied in [2] where in particular it is shown that the bifurcation structure for the stationary 
problem associated with (11.21) depends on the parameter e as well as on the length L of the 
interval Q. In fact, one can prove 

Proposition 1.1. JH, Proposition 3.1] The k-th bifurcation from the trivial solution of the sta- 
tionary problem associated with lil.^) is a supercritical pitchfork if L > kir/y^ and a subcritical 
pitchfork if the inequality is reversed. 

Through an analysis of the time map (see [10]) associated with the stationary problem for (II. 2p . 
one can also show that, for any given value of L, there is a value of e = €*{L) such that for 
e < e*{L) there are classical (i.e. C^((0,L)) nC^([0,L])) solutions to that problem and that 
solutions at e = €*{L) develop infinite gradient (see [2] for details). 

Since the physical arguments of Rosenau [7] for assuming the free energy of a system to have the 
form (11.11) are compelling, it makes sense to study a mass- conserving version of (11.21) . There are 
many ways of doing that, the simplest being to consider the mass-constrained L^-gradient flow 
of (II. ip which leads to 

ut = e{^{u^))^ + f{u) -7^ f{u) dx, xen, t>0 (1.3) 

l"l Jn 

where \Q\ is the length of Q and e, ijj{s) and f{u) are as above, with the Neumann boundary 
conditions = on dQ and a suitable initial condition u{x, 0) = Uq{x). Clearly, 



Thus (II. 3p is the quasilinear analogue of the Rubinstein-Sternberg equation from [8] . Local mass- 
conserving versions of (II. 2p can also be introduced leading to a Cahn-Hilliard type equation. 



ip{s) —7- ±1 as s — 7- ±oo. 




ut = -{e{tp{u^))^ + /(n))^.^., X eQ, t>0, 
Ux = {e{ij{ux))x + f{u))x = 0, X e dVt, 



(1.4) 
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which can be obtained as the iy"^-gradient flow of the free energy functional in (11.11) . Note that 
one can integrate the stationary problem associated with (11. 4p twice with respect to x using the 
Neumann boundary conditions to obtain the stationary problem associated with (II. Sp . Hence 
stationary solutions to (II. Sp are also stationary solutions to (II. 4p . 



2 The stationary problem 

We are interested in this section in characterising the multiplicity of solutions to the stationary 
problem associated with (II. Sp and hence also with (JL 

A '•^ 



mUx)), + Xf{u)--I f{u{x))dx = 0, XG(0,L), (2.1) 







1 '•^ 



Ux = 0, atx = 0, ,L, — / u{x)dx = M, 







as the bifurcation parameter A = - varies in (0, oo) 



Liapunov- Schmidt reduction 

We use the Liapunov-Schmidt reduction [1] to obtain local bifurcation results for (classical) 
solutions to the non-local equation in (12. ip . This work is a quasilinear analogue to work done 
in [3] which considers the case where ipis) = s in (12. ip and uses local bifurcation and path- 
following methods to examine the changes in bifurcation diagrams of stationary solutions to the 
Cahn-Hilliard model of phase separation as the mass constraint is varied. 



Note that if u{x) is a solution to (12. ip then u{L — x) is also a solution to (12. ip and so bifurcations 
from the trivial solution u{x) = M of (12. ID arise as pitchforks. 

We set V = u — M and recast problem (12. ip as G{v, A, M) = 0, where 

G{v,X,M)=(^=^^==] +Xf(v + M)-j r f{v{x) + M)dx. (2.2) 

Hence we regard G as an operator G : D[G) d H ^ H where D{G) is given by 
D(G') = jt; G C2((0,L)) : t;'(0) = t;'(L) = 0, j i;(x)dx = o|, 



and H is the space 



H=\^we C((0, L)) : ^ w{x) c/x = o| . 
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The linearisation about the trivial solution = is given by 
(c?G)o,A,M-w = ^G{0 + hw,X,M)\h=o 



w,, + Xf'{M)w - - / f\M)w{x) dx 
^ Jo 

w,, + Xf'{M)w, 



since w G D{G). Hence ker((iG)o,A,Af is one-dimensional when A = = ^2j/|jv/) ^^'^ spanned 
by Vk = cos (^)- Thus in a neighbourhood of a bifurcation point (A^, 0) of (12. ip . we show that 
solutions of G{v, A, M) = on if are in one-to-one correspondence with solutions of the reduced 
equation h{\,y) = 0, ?/ G M, obtained through a Liapunov-Schmidt reduction. 

For the details of the computations, the reader is referred to the Appendix. For the derivatives 
of the bifurcation function h{X,y) evaluated at A = A^, y = we obtain from (15. 2p 

h = hy = hyy = JlX = 0, 

Also, from ( 15. 9p . hxy = Lf\M) and hyyy is given by (15.61) . 

Let us consider first the case M = 0. Then h\y = L > by (15.90 and from (15. 6p 

Kyy = ^^(3A;V -6L2), 

so that the bifurcation from the trivial solution will be subcritical if L < ^ and supercritical if 
L > This corresponds with the result mentioned in Section [T] and obtained in [21 Proposition 
3.1] for the non-conserving Neumann stationary problem for (II. 2p . 

Now suppose < |M| < Using (15. 6p . we see that hyyy > as long as 

L<L*:= — \ \ 2.3 

which has a vertical asymptote when \M\ = and a turning point when |M| = We plot 
the relationship between L* in (12. 3 P and M for some G N in Figure 12.11 
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Figure 2.1: Plot of length L* against mass M for some A; G N, {k = 1). 



Also, we note that for all L > with -i= < |M| < 

3en'[f'{M)f + L'f"{M)f{M) + ^[f"{M)f 
= 3A;27r2(l -3M2)2 + 6L2(5M2 - 1) > 0, 

so that from fl5.6p . hyyy > for all for all such L and M. Therefore, since h\y = f'{M)L > for 
all M such that \M\ < we have established: 

Proposition 2.1. For < \M\ < the k-th bifurcation from the trivial solution of l{2.1\) is a 

supercritical pitchfork if L > L* and a subcritical pitchfork if L < L* (with L* given in Ii2.3\) ). 
For < \M\ < bifurcation from u{x) = M is a subcritical pitchfork for all L. In the 

parameter regime \M\ > there are no bifurcations from the constant solution u{x) = M for 
any L. 



Remark: Note that, unlike in the semilinear case of [5] wherein ip{s) = s, one can have both 
supercritical and subcritical pitchforks for different values of k. 
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3 Nonexistence of classical solutions for A large enough 

In this section we show that as A — )■ oo, branches of classical stationary solutions do not persist 
indefinitely. To do that, we introduce an auxiliary problem. For fixed L and a G M, consider 



u 



u{0) = u\L) = 



A/(m) = Xa, (3.1) 



which can be written as the first order system 



u = V 

v' = A(a-/(n))(l + t;2)f 



2.. , (3.2) 



with the first integral given by 



Hiu.v) = + Xiau - Fiu)). (3.3) 

V 1 + f 

The connection between (12. ip and (13.11) is seen by integrating (13. ip over Q = {0,L) to obtain 

/(n(x)) dx = a. 



1 '•^ 



L 







Lemma 3.1. For a fixed L and any a G M such that \a\ < there is a number A*(a, L) such 
thatWX > X*{a,L), there are no non-constant classical solutions to ( I A*. 



Proof. Suppose < a < so that the equation f{u) = a has three solutions ui, c, Ur with 
ui < < c < Ur where (m/, 0) and (m^, 0) are saddle points for (13.21) and (c, 0) is a centre. From 
considerations on the first integral (13. 3 p associated with the ancillary problem (13. ip one can see 
that for A < A/i(a), there is a homoclinic loop connecting {ur,0) to itself which, in the case of 
a particular A < A/i(a), we have denoted by 7a in Figure [XT] Note that for each a, Xh{a) is 
obtained through solving 

H{ur, 0) = H{c, —oo), 

for A = Xh{a). Non-constant classical solutions to (13. ip for A < A/i(a) are represented in the phase 
plane by trajectories which encircle the centre (c, 0), start and end on the M-axis and which are 
contained within the homoclinic loop 7a. 
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Figure 3.1: Phase portrait for with A < Xh{a) (a = 0.1, A = 3 < A/,(0.1) = 6.3426). 

For A large enough, i.e. for A = Xh{a), the homochnic loop will break. Through further consid- 
erations on H{u,v), one sees that if A > A/i(a), as in Figure l3T2| there exist values u*{X) and 
M**(A) such that u*{X) < c < u**{X), 

/7K(A),0) = if(w**(A),0), (3.4) 

and 

M*(A) c", M**(A) c+ as A ^ oo. (3.5) 
Note that for each A > A/j(a), we can find u*{X) by solving 

H{u*{X),0) = H{c,oo), 

for ui < u*{X) < c and then obtain u**{X) via (13. 4p . Hence non-constant classical solutions 
to (13. ip are now confined to 7^ as in Figure 13.21 where 7^ is the region in the phase plane enclosed 
by the trajectories through (m*(A),0) and through (m**(A),0) for A > Xh{a). 
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Figure 3.2: Phase portrait for (E^) with A > A/,(a), (a = 0.1, A = 15 > Aft(O.l) ~ 6.3426). 



Now assume that for a given A and a given n E N, there exists a non-constant classical solution 
u{x) to fl3.ip which has n points of inflection in (0,L). Such a solution will have period equal 
to — . For A large enough, because of (13. 5p . the only possible non-constant Neumann classical 
solutions to (13. ip will be in a small neighbourhood of the centre, hence we linearise (13. ip around 
(c, 0) to obtain 

u" + Xf'{c)u = 0. 

Thus the period of such a solution is equal to 

27r 



(3.6) 



which is well-defined since < c < ^ so that /'(c) > but (13. 6p will not be equal to ^ if A 
is large enough. This contradiction proves Lemma 13.11 in the case that < a < The case 
— < a < can be treated similarly. Note that in the case that a = (for which non-constant 
classical solutions to (13. ip are also stationary solutions of (II. 2p ). there are heteroclinic loops 
connecting saddle points which break for A large enough rather than homoclinic loops but the 
arguments carry through just as easily. □ 



Remark: We note that for each < a < if A is small enough then the phase portrait 
for (13. 2p contains non-trivial classical solutions to the ancillary problem (13. ip which are also 
non-trivial classical solutions to (12. ip for some < M < 1. Similarly, for each —-^^ < a < 0, if 
A is small enough then the phase portrait for (13.21) contains non-trivial classical solutions to (13. ip 
which are also non-trivial classical solutions to (12. ip for some — 1 < M < 0. If a = then for A 
small enough, the phase portrait for (13. 2p contains non-trivial classical solutions to (13. ip which 
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are also non-trivial classical solutions to ( 12. ip with M = 0. 



We now concentrate on the multiplicity of monotone classical solutions to (I2.ip where without 

loss of generality, < M < so that from now on, for a G 0, X*{a,L) represents the 

value of A for which there are no monotone classical solutions to (13. ip VA > A* (a, L). Note that 
A*(a, L) — A*(0, L) as a — and since the value A*(a, L) must occur after the homoclinic orbit 
has broken for a fixed a, we also have that A* (a, L) — )■ oo as a — 



Let us drop the dependence of A*(a,L) from Lemma [3.11 on L and denote the value simply by 
A*(a) since we will regard L as being fixed. Hence Lemma IXT] implies that for every a G 0, j ? 

there is a value A* (a) such that monotone solutions to the ancillary problem (13. ip develop infinite 
gradient at some xq G (0, L). Therefore for every A > A*(0), there is a value of a, corresponding 
to the inverse of A* (a), which will be denoted by a* (A) and be such that a* (A) — > as A — j- oo 
and a*(A) — )■ as A — )■ A*(0)"''. We give a sketch of the function a*(A) in Figure 13731 but note that 
we do not prove that the curve a* (A) is continuous (see Theorem I3.3p . 




Figure 3.3: Sketch of the curve a* (A) in the (A,a)-plane. 

Lemma 3.2. For a fixed L and each < M < the bifurcation curve of monotone classical 
solutions to l\2. 1\) does not exist in the region in the {X,a)-plane defined by a*(A) < a < for 
A large, where a* (A) is the curve obtained from Lemma \3. 1\ 



Proof. Let a(A) denote the bifurcation curve of monotone classical solutions to (12. ip in the (A, a)- 
plane and assume the contrary to Lemma [3.2[ Then there exists a sequence (A„)^]^ such that 
A.„ — > oo as n — > oo and 

a*(A„) < a(A„) < ^ 



3v^' 

for all n G N sufficiently large. Then, since a*(A„) — )■ as n — oo, we have by the sandwich 



Quasilinear bistable equation 



10 



theorem that 

2 

fl(A„) — 7- — -1= as n — 7- oo. (3.7) 
3v 3 

But M is constant along a(A) and (13 .Zp imphes that M must be equal to which is a contra- 
diction. □ 

We now have the following theorem regarding (monotone) classical solutions to the non-local 
mass-conserving equation (12. ip : 

Theorem 3.3. For fixed L and < M < there exists a value Xi{M,L) such that for 
A > Ai(M, L) there cannot exist monotone classical solutions to the non-local mass- conserving 
equation in h2.1\) . 

Proof. In the (A,a)-plane, the curve a* (A) obtained from Lemma l3.ll separates two regions; a 
region in which monotone classical solutions to the non-local mass-conserving equation (12. ip can 
exist for fixed L > and < M < and a region in which monotone classical solutions 
to (12. ip cannot exist for such L and M. As in Lemma [3. 2 [ let a(A) denote the bifurcation curve 
of monotone classical solutions to (12. ip . We know by Lemma [3.21 that a(A) must intersect the 
curve a* (A) at some point since a(A) cannot not exist in the region defined by a* (A) < a < -^j^ 
in the (A, a)-plane for sufficiently large A and from the Rabinowitz theorem (see [HI Theorem 
13.10]), it has to go somewhere as A increases. We want to prove that a(A) intersects a*(A) at a 
point of continuity of a* (A). 




Figure 3.4: The curve a* (A) and a proposed bifurcation curve a(A) of classical monotone solutions 
to (12. ip intersecting a* (A) at a supposed point of discontinuity of a* (A). 



We assume the contrary and consider the bifurcation curve a(A) of monotone classical solutions 
to (12. ip in the (A,a)-plane as we have it in Figure 13.41 We assume that the curve a* (A) is 
discontinuous at some value A^ and that the bifurcation curve a(A) enters the region of no 
classical solutions to (12. ip at this point of discontinuity of a* (A). By the Rabinowitz theorem [9l 
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Theorem 13.10] there must be a classical solution in a neighbourhood of this point and the 
bifurcation curve of classical solutions can be extended. However, we would then have entered 
into the region in which there can be no classical solutions to (12. ip and so we have obtained a 
contradiction. Hence the bifurcation curve must intersect the curve a* (A) at a point of continuity 
and the theorem is proven. 

□ 



Remark: Although we established Theorem 13.31 for monotone classical solutions to (12. ip . the 
result can be generalised to non-monotone classical solutions u{x) to (12. ip which have n points 
of inflection in (0, L) for some n G N. Thus for each n G N, there exists a value A„,(M, L) such 
that for A > A„(M, L), there are no classical solutions to (12. ip with n points of inflection in 
(0,L). However, we point out that unlike the situation in the (semilinear) Rubinstein- Sternberg 
equation [3] , there is no obvious way to deduce the behaviour of a particular branch of solutions 
to (12. ip from that of the monotone one. Also, in the proof of Theorem 13.31 we did not establish 
that the curve a* (A) is continuous but we did show that it must be continuous at the point A^ 
where the bifurcation curve of monotone classical solutions to (12. ip intersects it. 



We have not given a precise depiction of the bifurcation curve of monotone classical solutions 
to the mass-conserving problem (12. ip in the (A,a)-plane for given L and < M < Hence 
for fixed L and M subcritical and then for fixed L and M supercritical, we numerically obtain 
the curve a* (A) from Lemma [3.11 and the bifurcation curve a(A) of monotone classical solutions 
to (12. ip and plot these in the (A,a)-plane. With regards to problem (12. ip . we know that along 
the line of trivial solutions, = 7; /g^ fi'^i^)) dx is constant and equal to M — M^. Also, given 
L, for each < M < bifurcation from the trivial solution of ( 12. ip occurs when A = ^2^(j^^-) 
which we can rearrange to 



1 

M 



and since bifurcation points appear along the line of trivial solutions (upon which a = M — M^) 
we can plot a curve of bifurcation points for fixed L in the (A, a)-plane by considering the function 

ab{X) = M-M^ 



1 71^ 



3v^V L^X \ L^X, 

For a given L, we can numerically plot the curve a* (A) in the (A,a)-plane by solving (12. ip for a 
variety of M G 0, ^ j using path-following methods of AUTO [I] and determine for each M the 

value of A and the value of a = -^^ f{u) dx for which the stationary solutions to (12. ip develop 
infinite gradient in (0, L) and then plot a against A. We can then use AUTO again for fixed L 

and a particular M G 0, to plot the bifurcation curve a(A) of monotone classical solutions 

to (12. ip in the (A, a)-plane and obtain the value Ai(M, L) of Theorem 13.31 
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Thus we fix L = 1 and take M = 0.1 so that the bifurcation at A = -pjQj^ = 10.1749 from the 
trivial solution is subcritical and a = 0.099 along the line of trivial solutions to (12. ip . According 
to AUTO, the branch of monotone classical solutions to (12. ip for these parameter values stops 
when A = Ai(0.1,l) = 5.6579 with 5(Ai(0.1,l)) = a*(Ai(0.1, 1)) = 0.0289 and we have plotted 
all relevant curves for these values of L and M in the (A, a)-plane in Figure 13.51 

a 



2 

3v^ 




Figure 3.5: The bifurcation curve a(A) of monotone classical solutions to (12. ip for L = 1, M = 0.1 
and the curves ab{X) and a*(A) in the (A,a)-plane. 

Now consider L = 2.5, M = 0.3 for which a = 0.273 along the line of trivial solutions to (12. ip . 
In Figure [331 we plot all relevant curves in the (A, a)-plane for these values of L and M. We have 
a supercritical bifurcation from the trivial solution when A = 2 5^/'(o 3) ~ 2.1632 and according 
to AUTO, the branch of monotone classical solutions in this case stops when A = Ai(0.3,2.5) = 
4.0860 with a(Ai(0.3, 2.5)) = a*(Ai(0.3, 2.5)) = 0.0051. 



Quasilinear bistable equation 



13 




Figure 3.6: The bifurcation curve a(A) of monotone classical solutions to (12. ip for L = 2.5, 
M = 0.3 and the curves afe(A) and a*(A) in the (A,a)-plane. 



On the basis of the preceding numerical experiments we would conjecture that the bifurcation 
curve a(A) of (monotone) classical solutions to (12. ip is monotonic and for given L and M, a is 
at its largest along the line of trivial solutions i.e. when a = M — M^. 

Remark: We have not discussed the case of having < \M\ < 1 in (12. ip . i.e. what happens 
in the "metastable" regime in which there are no local bifurcations from the trivial solution 
and for which the trivial solution is always stable. As \M\ — )■ f'{M) — )■ and the first 

and all subsequent bifurcation points A^ = ^.^/'Im) S° infinity. In the semilinear situation 

studied in [3], in passing from the spinodal to the metastable regime, they show through spectral 
approximations and path-following methods that the saddle-nodes which exist for < \M\ < 
move off to the right as f'{M) — )■ O"*" but at a speed much slower than that of the bifurcation 
points. In our case, for all L with M just less than ^ so that the bifurcation from the trivial 
solution is subcritical by Proposition 12. H the classical solutions to (12. ip stop existing before we 
reach a saddle-node. Hence we were not able to perform a two parameter continuation in A and 
in M of the saddle- nodes for \M\ beyond We can however say what happens in the parameter 
regime \M\ > 1. By the remark after Lemma [3.11 we see that it is not possible to construct a 
non-trivial classical solution to the ancillary problem (13. ip which will have average mass M such 
that \M\ > 1. One can also see from phase portraits associated with (13. 2p that it is also not 
possible to construct a non-classical solution to (13. ip for any \a\ < which will have \M\ > 1. 
Therefore for \M\ > 1, there are no non-trivial solutions to (12. ip for any A G (0, oo) which is also 
true of the semilinear problem (see [3]). 
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4 Numerics 



In this section we carry out some numerical experiments for the non-local mass-conserving equa- 
tion in fll.3p in order top see what (stable) non-classical solutions of fll.3p look like. We first 
derive an explicit numerical method to solve the equation (11. 3p and we are grateful to Dr John 
Mackenzie of The University of Strathclyde for discussions on the numerical scheme we outline 
in the following subsection. 



Numerical approximation 

We derive a mass- conserving numerical scheme to solve 



Ut 



l+«2 



X 



Ux{0,t) = Ux{L,t) = 0, {x,t) e dn X {0,T) (4.1) 
u{x,0) = Uq{x), X E fl, 

where X e {0,oo), Q = {0, L), L > 0, f{s) = s - and < T < oo. We obtain (gl]) from fOD 
by multiplying (II. Sp by ^ and scaling time as t h-> et. One can easily show that the non-local 
equation (14. ip conserves mass and to hope to have the same at the discrete level we must discre- 
tise the equation in conservative form, i.e. in the form given in (14. ip . 



We discretise the space interval into + 1 evenly spaced points 

= Xi < a;2 < . . . < xn+1 = L, 

so that Ax = jj^ and we regard there as being cells of width Ax around each internal 

point Xi, i = 2, ■ ■ ■ , N while at the boundaries i = 1 and z = -|- 1 we use cells [xi, xa] and 
[a^AT+i, xn+i] of half width ^ as in Figure ITTl 

L 



xix^_iX2 xs X4 ^N^iXjM Xj^^XN+i 



Figure 4.1: Discretisation of the space interval [0,L]. 



At a particular Xi E fl = (0, L) (hence i E {2 ■ ■ ■ N}), tn G (0, T) the equation is given by 

Ut{Xi,tn) = ( —^^^====] + Xf{u{Xi,tn)) - Y [ f{u{x,tn))dx. 
\y/l + ul{Xi,tn) J ^ LJo 
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Let v{x,t) denote the flux, i.e. 



v{x, t) 



so that at some Xi E fl = (0, L), tn E (0, T) 



1 ~1~ li.j^iyXi^ tji^ 



Ax 



(4.2) 



that is, we approximate the flux term by a first central difference in space at the point Xj with 
half-spacing. 



Now, 



v{Xi^l,tn) 



l + ul{Xi^l,tn) 



{u{Xi+i, tn) - u{Xi, t„))/Ax 



u{Xi-^.l,t„)-u{xi,tn) 

Ax 



(4.3) 



and 



v(Xi_i,tn) 



UxiXj 1 , tn) 



1 + ulix. 



{u{Xi, tn) - u{Xi^i, tn))/Ax 



u{Xi,t„)-u{Xi^l,tn) 

Ax 



(4.4) 



so that in (14.31) we approximate the term f (Xj_,_i , t„) by a first central difference in space about the 
point with half-spacing and in (14. 4p we approximate the term f(Xj_i,t„) by a first central 
difference in space about the point x^_i with half-spacing. Therefore by (14. 2p . (14. 3 p and (14. 4p 



^x (-^ii tn) 



V^l + ul{Xi,tn) 



'^x i,Xi , tn) 
( 



1 



Ax2 



(^(Xi+i, tn) - u{Xi, tn)) {u{Xi, tn) - n(Xi_i, t„)) 



u{Xi-f-i,tn) — u{Xi,tn) 

Ax 



u{Xi,t„)-u{Xi^l,tn) 

Ax 



\ J 



Quasilinear bistable equation 



16 



We approximate the integral term in (14. ip at some given time t„ using the midpoint rule for 
numerical integration as follows 



/(n(x,t„)) dx 



^3 
1 



f{u{x,tn))dx+ f{u{x,tn))dx+ ...+ f{u{x,tn))dx 



1+^ 



+ ...+ / f{u{x,tn))dx+ / f{u{x,tn))dx 

Ax „, , v-^ . „, , Ax 



'-f{u{xi,t^)) + ^ Ax/(u(xj,t„)) + — /(n(xAr+i,t„)) 



i=2 



Let u{xi,tn) = so that using the forward Euler method, we approximate (14.11) at interior 
points Xj G (0, L), « = 2, ■ ■ ■ , by the numerical scheme 



+At(^ 



1," — 



1/", , -u" 



z z — 1 



+ A/«) 



(4.5) 



for i = 2, . . . , N. At the boundary points x = and x = L we consider the discretisation in the 
half-cells [xi,X3] and [x^_,_i , xat+i] respectively, so that we have 



and 



At 



A 
L 



Un — 



Ax' 



1 + 



[ Ax 



N 



i=2 



(4.6) 



7/*^+^ —7/" 



At I 



A 
L 



u 



Ax2 



,n _,,n 



Ax 



N 



i=2 



(4.7) 



For the above numerical scheme to approximate (14. ip reasonably it must also conserve mass. 
Therefore we have the following theorem. 



Theorem 4.1. The explicit numerical scheme contained in \4-5\j, \4-()\j, ([^. 7[) conserves mass. 
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Proof. We multiply (14 .Sp by Ax and sum over i = 2,3, ■ ■ ■ ,N to obtain 



i=2 i=2 \ 1=2 



K+1 - < 



1 + 



,71 \ 2 



Ax 



Ax 



N 



N 



Xj2fK)Ax-J2Aa 



i=2 



i=2 



L 



N 



i=2 



N 



y u"Ax + At — 
^ I Ax 



i=2 



TV 



Ax 



1 + 



1 2 



Ax 



a5^/«)Ax-EAx 



i=2 



A 



i=2 



(4. 



Now we multiply both of the boundary terms in (14.6^ and (14. 7p by and add the resulting 
equations to f l4.8p to give 



u 



n+l 



N 



i=2 



+1 I "-N+l 



Ax 



i=2 



Ax 



and so mass is conserved at the discrete level as required. 



□ 



Numerical experiments 

We present the results of some numerical simulations for the non-local mass-conserving equa- 
tion (14. ip using the explicit mass-conserving numerical scheme described in Section HI For 
< M < fixed, we choose both sub- and supercritical lengths L (see Proposition 12. ip . 
In the subcritical case (Experiment 14. 2p we use AUTO to plot the bifurcation diagram for mono- 
tone classical solutions to (12. ip and find the value Xi{M,L) of Theorem 13.31 for these values of 
M and L. Then we solve (gl]) using (03]), (gj]), (HID for A < Xi{M,L) and for A > Xi{M,L) 
with monotone initial data satisfying the mass constraint and present the initial data with the 
final equilibrium state in each case. In the supercritical case (Experiment 14.30 . we use AUTO to 
plot the bifurcation diagram for monotone classical solutions to (12.11) and also for non-monotone 
classical solutions to (12. ip which have two inflection points in (0,L). We find Xi{M,L) and 
A2 (M,L) the values of A such that for all A > Xi{M,L), there are no classical solutions to (12. ip 
with i inflection points in (0, L) for 2 = 1,2 respectively. We then run experiments solving (14. ip in 
the cases A > Ai(M, L) and A > A2(M, L) using (14. 5p . (14. 6p . (14. 7p for particular initial data uo{x). 
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Experiment 4.2. We fix M = 0.2 so that L* = 2.1856 from (12. 3p and so, to have a subcritical 
L, we can take L = 1.7 in this case. We use AUTO to plot the bifurcation diagram for monotone 
classical solutions to (12. ip for these values of L and M in Figure 14.21 We have a subcritical 

2 

bifurcation from the trivial solution u{x) = M when A = i 72 /'(o 2) ~ 3.8808 and according to 
AUTO, the bifurcation diagram for monotone classical solutions to the stationary problem stops 
when A = Ai(0.2, 1.7) ~ 4.3032 as in Figured 

«(0) 



AI - 
0. - 




Figure 4.2: Bifurcation diagram for monotone classical solutions to (12. ip with M = 0.2, L = 1.7 
where Ai (M,L) ~ 4.3032. 



We take the following initial data 

uo{x) = 0.3 - 0.5 tanh (1000 
which satisfies the mass constraint, i.e. 



1-0,4 



(4.9) 



y [ Uo{x) dx = M = 0.2. 
L Jo 



Hence we present the equilibrium solutions to the time- dependent problem (14. ID . (14. 9 p in Fig- 
ure |4]3] for A = 4 < Ai(M, L) (left) and A = 5 > Ai (M,L) (right) obtained using the scheme 
in (iS!), (iJD, (iTD with = 500. 
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X < (0.2, 1.7) 



X' 




Figure 4.3: Initial data and equilibrium solutions to (14. ip . (14. 9 p for A < Ai(0.2, 1.7) (left) and 
A > Ai(0.2,1.7) (right). 



Experiment 4.3. Suppose we now consider the case M = 0.2 with L 
supercritical bifurcation from the trivial solution u{x) = M when A 



2.5. We have a 
2.5y'{o.2) = l-'^945 and 

according to AUTO, the bifurcation diagram for monotone classical solutions to the stationary 
problem stops when A = Ai(0.2,2.5) ~ 4.0433. There is also a subcritical bifurcation from the 
trivial solution when A = ^ ^^Jrf^Q 2) ~ 7.1779 and a curve of non- monotone classical solutions 
to (12. ip which stops when A = A2(0.2,2.5) ^ 4.9872 just after it has reached a saddle-node at 
A = A^„ ~ 4.9714 as in Figure |4]4l 



uiO) 



M 

0.- 




10. 



Figure 4.4: Bifurcation diagram for classical solutions to (12. ip with at most two inflection points 
in (0,L) where M = 0.2, L = 2.5, Ai(M,L) ~ 4.0433, A2(M,L) ~ 4.9872 and there is a saddle- 
node at A = A,„, ~ 4.9714. 
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Suppose we solve (14.11) for L = 2.5, M = 0.2 and fixed A > Xi{M, L) with initial data given by 



uo{x) = /3 + 0.5 tanhl 1000 



7 



(4.10) 



and we vary the position xq = 7L of the interface in (I4.10p but ensure that 



1 
L 



Uo[X) 



(4.11) 



still holds by changing /3 accordingly. The results of taking L = 2.5, M = 0.2 with A = 8 > 
Ai(0.2,2.5) ^ 4.0433 and solving (14. ip for various values of 7 (and (3) are plotted in Figure H75] 
where one sees that modifying the initial data while ensuring that (14. lip holds has an effect 
on the equilibrium state to which the solution converges as t — cxo which is not the case for 
A < Ai(M, L). This suggests that the discontinuous equilibria for (II. 3p existing for A > Ai(M, L) 
are normally stable in the sense of [6]. 



0.5 



-0.5 



-y=0.5, |3 = 0.2 
7=0.6, |3 = 0.1 
-y=0.9, p = -0.2 



1.25 

X 



2.5 



Figure 4.5: Equilibrium solutions to (jH]), flCTj) with (7,/?) = (0.5,0.2), (0.6,0.1), (0.9,-0.2) 
solved using (jts]), dM]), (gT]) with N = 500. 



Finally we take A > A2(0.2, 2.5) with the following non-monotone initial data we define piecewise 
as follows 



Uo{x} 



f -0.2] ) 
^ -0.8 " 



) 





L 
2 



<X<f 

< X < L 



(4.12) 



0.32- 0.6 tanh (1000 
0.32 + 0.6 tanh (1000 

which satisfies the mass constraint (14. lip . In Figure I4l6] we show this initial data and the final 
equilibrium state to which the solution converges as t 00 for A = 8 > A2(0.2, 2.5). 
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Figure 4.6: Non-monotone non-classical equilibrium solution to f l4.ip . fl4.12p for A > A2(0.2, 2.5) ~ 
4.9872 solved using g3]), gS]), gZD with = 500. 



5 Conclusions 

We have introduced a non-local mass- conserving quasilinear equation (11. 3p that can be viewed 
as a one-dimensional quasilinear version of the Rubinstein-Sternberg equation from [8] and we 
have considered its associated stationary problem. We have shown in Proposition 12.11 that un- 
like the situation for the bifurcation diagrams of stationary solutions to the one-dimensional 
Cahn-Hilliard equation studied in |3] , the bifurcation behaviour for classical stationary solutions 
to (11. 3p as A is varied depends on the length L of the space interval Q as well as on the av- 
erage mass M of a solution. We have also proved (Theorem 13.31) that for each n G N there is 
a value of A = A„(M, L) such that for any A > A„(M, L) there cannot exist classical station- 
ary solutions to (11.31) with n points of inflection in (0,L). We have presented some numerical 
results using AUTO which illustrate Proposition 12.11 and Theorem 13.31 and we derived a mass- 
conserving numerical scheme for (II. 3p . There is numerical evidence (see Figure S3]) to suggest 
that for A > Xi{M,L), there exist a continuum of normally stable (in the sense of [6]) non- 
classical monotone stationary solutions to (11.31) with a discontinuity at some Xq G (0,L). We 
have also presented numerical evidence which suggests that there are also stable non-monotone 
non-classical stationary solutions to (I1.3P for A large enough (see Figure l¥3|l . 

An understanding of exactly what happens in the "metastable" regime (^-^ < \M\ < 1^ is re- 
quired. We were not able to perform a two parameter continuation in A and in M of the 
saddle-nodes arising in the < 1^1 < 75 parameter regime since for M large enough, the 
classical solutions to ( 12. ip stop existing in the bifurcation diagrams before the saddle-nodes are 
reached. One can show however through an analysis of the phase plane associated with the 
ancillary problem for (12. ip that for \M\ > 1, there are no non-trivial solutions to (12. ip . 



We also note that no existence or stabilisation theory for (II. 3p or (II. 4p is available. 
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We hope that the class of models introduced in this paper will be of interest to material scientists 
and that its applicability to real materials will be tested as has been done for more classical models 
in, for example, [5]. 



Appendix 



In this appendix, we present the details of the Liapunov-Schmidt Reduction we discussed in 
Section [3 For the operator G : D{G) H given in let S = {dG)o,x^^M ■ D{G) H, with 

kernel K. and range TZ given by 



span 



cos 



kirx 
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and TZ 



w e H : 



w{x) cos 



kfix 



dx = 



(5.1) 



respectively, and let E : H TZ denote the projection of H onto TZ. The linearisation operator 
S for equation (12.11) at a bifurcation point (Afc, 0) is a self-adjoint elliptic operator and therefore 
Fredholm of index zero (see [H Appendix 4]) and so, following the steps in [H p. 29], the spaces 
D{G) and H are decomposed as 



since /C 
then 



and /C^ 



D(G)=/C©/C^, if = /C©/C^, 

TZ. The coordinates chosen in the Liapunov-Schmidt reduction are 



2 cos 



/ knx 



Vk 



cos 



/ kTix\ 



From [U Chapter 1], the derivatives of the reduced function h{X,y) evaluated at A = A^, y = 
are given by 

hy = {vl,dG{vk)) = 
Ky = {v*kyd^G{vk,Vk)) 

hyyy = {vld'Givk,v,,Vk)-3d^Givk,S''E[d^Givk,Vk)])) (5.2) 
hx = {vl,Gx) 

hxy = {vl,dGx{vk) - d^Givk,S-'EGx)) , 
where (■, ■) denotes the L^-inner product on [0,L]. 



We now proceed to calculate the derivatives of the reduced function h{X, y) in order to determine 
locally the direction of the pitchfork bifurcations from the trivial solution u{x) = M of (12. ip for 
a given L and a given M. This is complicated by the fact that the operator G in (12.21) is not 
odd and so terms involving 5*"^ in (15. 2p will not necessarily vanish. In order to invert 5* we will 
need to solve an ordinary differential equation; see (15.30 . 



Quasilinear bistable equation 



23 



We have 



52 

((i^G')o,Afe,M(wi, W2) = Tr-^G{0 + tiwi + t2W2, Afc, M) 



dtidt 



L Jo 



and so 



Xkf"{M)w,W2 - / r(M)ii;i(x)w;2(x) rfx, 



''r(M)cos2 (^] dx 



Xkf'iM) cos' 



\ L J L J, 



L 



Hence by 



hyy = {vl,d^G{vk,Vk)) 
0. 



X 2 I kirxX ^ f"{M) 



da; 



Set 



so that 



and 



{d'G)o,xM^k,v,) = X,nM) cos' ( ^] - \,tiMl=p{x). 



p\0)=p\L) = 0, 



— / p{x) dx = 0, 



therefore (rf^G)o,Afc,Af (^^fc, ^^fc) e H. 



However, 



{d G)o^Xt:,M {vk,Vk) COS 



Xkf'iM) cos' 



[kTix\ f"{M) [knx 



X, 



■ cos 



dx 



so that {d''G)o,\f,,M{vk,Vk) G 7?. by (15. ip . Therefore {d''G)o,Xf,,M{vk,Vk) G trivially decomposes 



as 



((i^G)o,Afe,A/K,ffc) = (c?^G)o,Afe,A/K,^^fc) +0, 
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where {(PG)o,\,.,M{vk,Vk) G 71 and of course G 71^. Thus, since E : H ^ 71 is the projection 
of H onto the range of S, we have 



and we consider 



S~^E[{d'^G)o^Xk^Mivk,Vk)] = S~^{d^G)Q^XkM{'"k,Vk) = l{x) 

=^ {d'^G)o,Xk,MiVk,Vk) = Sl{x). 

Thus the second order ordinary differential equation that we need to solve for l{x) in order to 
obtain a formula for S~^E[{d'^G)o^Xk,M{vk,Vk)] is 



l"{x) + Xkf'{M)l{x) = \kf"{M) cos^ 
which has solution 

/(x) = cos 

so that 



(5.3) 



fknx\ lf"(M) f2k'nx\ 



ol^r. ,2^N . NT fkTlx\ lf"{M) f2kTTx\ 

S~'E[{d'G)o,xMvk,Vk)] = COS i^-j-j - Q^J^cos [-J^j 



Hence we compute 

d^G{v,,S ^'E[d^G{v,,Vk)]) 

/ knx 



cos 



d^G j^cos 



knx\ lf"(M) f2k'Kx\\ 

—)-inM-r\—)) 

fkTTx\ lf"{M) /2kTTx\ fkTix\ 



\ L ) 6 /'(M) V ^ 



cos 



1/"(M) /2fc7rx\ (k^x 



klTX 



L J Q f'{M) \ L J 



cos 



COS 



COS 



kTix\ lf"{M) (2kTxx\ fknx 



cos 



L J Q f'{M) \ L 



cos 



L 



L 



dx 
f"{M) 



and so we have 



{vl,3d^G{vk,S -'E[d^G{vk,Vkm 



6\kf"{M) 



cos 



k7rx\ lf"(M) {2kTix\ 

—)--,-TW)'"^\—) 



cos 



knx 



dx 



3\kf"{M) cos 



knx 



dx 



4 /'(M) 



AL [f'{M)f 



(5.4) 
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In addition to this, 



^3 

(d^G)o,Afc,A/(wi, W2, W3) = r., o, r,. G{0 + tiWi + t2W2 + t3W3,\k,M)\t^=t2=t3=0 

= -3[W'^W[W2 + W2w[w'^ + W'(W2W'3] + \kf"\M)wiW2Wz 
1 



-Xkf'"{M)— / wi(x)m;2(x)w3(x) rfx, 



1 '•^ 



so that 

(o!'G)o,A„MK,^;fc,^;fc) - 9[v'l{v',)'] + Xkf"iM)vl- XkriM)j- vi{x)dx 

- cos — — sm^ — — + Xkf (M) cos^ 



and 



{vl,d^G{vk,Vk,Vk)) 



Therefore from (15. 2p . (15. 4 p and (15. 5p we obtain 

^ 3A;27r2[/'(M)]2 + L^f"'{M)f{M) + — [/"(M)]^ 



Also, 

Gxiv,X,M) = fiv + M)-]r [ f{v{x) + M)dx, 

^ Jo 

so that Gx{0, Xk, M) = which imphes that 

(o!'G')o,A„M(t^fc, S"1E[Ga(0, Xk, M)]) = 0, 

while 

{dGx)o,x„M-y^ = f'{M)w-]r [ f'{M)w{x)dx 
= f'{M)w, 

for any w G D{G). Therefore, from (15.20 . (15.70 and (15. 8p we have 

hxy = (vldGxivk) - d'G{vk,S~'EGx)) 

= J^'2cos(^)/'(M)cos(^)rfx 
= Lf{M). 
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